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^ . Abstract 

vn ■ 

f~^ I Non-miniinally coupling a scalar field to gravity introduces an additional curva- 

ture term into the action which can change the general behavior in strong curvature 
regimes, in particular close to classical singularities. While one can conformally 
^^ I transform any non-minimal model to a minimally coupled one, that transformation 

^P ■ can itself become singular. It is thus not guaranteed that all qualitative properties 

O . are shared by minimal and non-minimal models. This paper addresses the classi- 

cal singularity issue in isotropic models and extends singularity removal in quantum 
gravity to non-minimal models. 
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X 1 Introduction 

It is well-established that general relativity implies space-time singularities for most solu- 
tions deemed to be relevant for physical situations such as cosmology or black holes [H |2] . 
Precise general properties of such singularities, besides the criterion of geodesic incom- 
pleteness usually used to prove singularity theorems, are however difficult to find. Even 
in isotropic models there are several types of singularities in addition to the common big 
bang type singularities at vanishing scale factor a{t) of a Friedmann-Robertson-Walker 
model 

ds' = -dt' + a{ty (-^^ + r^dnA . (1) 



Other possibilities include big rip singularities where the scale factor diverges in finite 
proper time, such as in superinfiationary solutions of constant equation of state parameter 
w < — 1, or sudden singularities at non-zero, finite a. These alternatives are rather special: 
While a big bang singularity can be shown to arise for isotropic models using only energy 
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conditions and non-zero expansion at one time, singularities at non-zero a require additional 
conditions to be satisfied by the matter ingredients. Such singularities can therefore be 
obtained in phenomenological matter models, but not so easily in models sourced by scalar 
or other fields. In minimally coupled scalar field models, for instance, one has only the 
freedom to choose a potential which does not allow, e.g., singularities associated with 
super-inflationary expansion. 

One additional freedom for scalar fields is a non-minimal coupling to space-time cur- 
vature by adding a term —^^^RcfP' to the Lagrangian. This can change in particular the 
behavior close to curvature singularities where R becomes large. Such a term is also in- 
teresting because it is the only possible curvature coupling of any matter field which does 
not require one to use additional length scales for the correct dimensions. Often, curvature 
couplings of fundamental fields are considered as effective terms to be derived from a quan- 
tum theory of gravity. In such a case, there is an additional length parameter provided by 
the Planck length £p = yGH which multiplies effective curvature terms. All these terms 
thus vanish in the classical limit except for the non-minimal coupling term of a scalar with 
quadratic coupling function, which is not multiplied by £p. Since a classical theory with 
such a term makes sense, one can also ask whether it would change any results of quantum 
gravity obtained with minimal coupling. Other curvature couplings, which are multiplied 
by powers of £p, cannot reasonably be quantized because they arise only from an effective 
description of an already quantized theory of gravity. 

One application of quantum gravity is singularity removal which can depend on the 
matter type. Moreover, with a non-minimal coupling term of a scalar already the classical 
behavior can change close to curvature singularities. We therefore perform in this paper 
an analysis of isotropic singularities in non-minimal scalar field models based on canonical 
transformations to a minimally coupled model. We will end with conclusions for singularity 
removal by quantum gravity. 

2 Non-minimal coupling 

Non-minimal coupling is usually introduced at the Lagrangian level which makes proper- 
ties of Lorentz invariance manifest. At the level of the action one can also see that any 
non-minimally coupled model allows a conformal transformation and re-definition of the 
field variables which brings it into minimally coupled form [3 II]. This transformation, 
however, can be difficult to determine analytically which makes the discussion of singu- 
larities complicated. After recalling these properties, we will thus switch to a canonical 
formulation. This allows us to perform canonical transformations which are more general 
than transformations allowed at the Lagrangian level. It will turn out that a much simpler 
transformation to minimally coupled form is then possible which we will exploit in our 
discussion of singularities. 



2.1 Equations of motion 

The simplest form of a non-minimally coupled scalar field action is given by 

S[gat, <P] = I d'x^/^ {^^R - ^(7"^V.0V,0 - f/(0)) (2) 

where g = det{gab), k = SnG, R is the scalar curvature of the metric and U{(f)) is the 
self-interacting potential of the matter scalar field 0. The curvature term is coupled to 
matter through the coupling function /(0) in a form which is manifestly Lorentz invariant. 
One can view this function as defining an effective gravitational constant k, in the sense 
f^es = f^/f{<P) and hence we require /(0) to be always positive. If / is thought of as 
arising from a perturbative expansion of a field theory, it is a polynomial and restricted 
by boundedness from below to be an even function. Finally, the coupling function should 
satisfy the limit /(0) -^ 1 when 0—^0, recovering the standard minimally coupled action. 
Given the restrictions above, the simplest non-trivial form of the coupling function is 

/(0) = 1 - a0l (3) 

In fact, this is the form of coupling function that is normally considered in the literature. 
We will refer to a = k^ as the coupling strength, vanishing a corresponding to the case of 
minimal coupling. We will later see, however, that the behavior around singularities can 
sensitively depend on the form of /(</>). 

Upon reducing the action to a flat Friedmann-Robert son- Walker model (^ it becomes 
a functional only of the scale factor a and the homogeneous scalar 0. Its canonical version 
depends on these two values and their momenta. In what follows, we will use independent 
variables p G M, such that \p\ = a^, and c = 7a. They have geometrical meaning as a triad 
component p (such that sgn(p) determines the orientation of the triad) and connection 
component c = 70, (or extrinsic curvature), which corresponds to using isotropic Ashtekar 
instead of metric variables El E! • The positive real number 7 is the so-called Barbero- 
Immirzi parameter j^l E] which we keep for generality but can be assumed to equal one 
here. These variables are the basis for a loop quantization to be employed in the final 
section. 

With these variables the action Q becomes 
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(4) 



with the lapse function A^. The dot denotes a proper time derivative. From (ji} we see 
that the pairs (c, /(0)p) and (0, vr = |p|'^/^A^~-'^0) are canonically conjugate variables with 
non- vanishing Poisson brackets 

{C,/(0)P} = ^, {0,77} = 1 (5) 

Note that the Poisson bracket {p, vr} = /(0)p{/(0)^^, tt} is not zero, in contrast to minimal 
coupling, while {/(0),p} is. This will motivate our choice of the new canonical variables 



transforming to a standard minimally coupled form later on. The dynamics is generated 
by the Hamiltonian H = N{Hg + H^f,) where 



It is now very convenient to proceed to conformal time by setting the lapse N := ^/\p\. 
The Hamiltonian constraint then becomes 

2.2 Conformal transformation 

A conformal transformation Qab = f{4>)gab transforms the action with non-minimal coupling 
to 



with 
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for a coupling function of the form (jS)). This can be brought to the standard minimally 
coupled form for the redefined scalar = / d0F(0), implying a rather complicated relation 
0(0), e.g., 

76^0+ Vl- ^(1-600' 




;i - 60 arcsin ( Va(l - 6^0 ) + \/ 2 ^°S 
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ife<l/6. 

From the canonical point of view, moreover, a conformal transformation implies a more 
complicated change in the momenta conjugate to the spatial metric. As we will see now, 
however, the larger freedom in canonical transformations compared to conformal ones 
allows a similar transformation to minimally coupled form which even leads to a simpler 
redefinition of the scalar. This will be exploited in the rest of this paper to discuss the 
behavior of isotropic singularities which will be much more transparent than discussions 
based on conformal transformations of the space-time metric. 

2.3 Canonical transformation 

In order to perform a canonical transformation bringing the Hamiltonian into its standard 
(i.e. minimally coupled) form, we start with the geometrical variables {c,p) and define 

g := c, p := /(0)p (8) 



such that {c,p} = ^. Note that this does not correspond to a conformal transformation 
from the original gravitational pair {p, c) to a new one {p, c) because in that case c would 
have to change, too. The Hamiltonian thus takes the form 

K72 2\p\ f[(t)f 

In this form, the Hamiltonian is written in terms of canonical variables (p, c, 0, tt) and 
would thus be the starting point for a canonical quantization. Finally, redefining the field 
variables 
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such that {</), vr} = 1 we obtain 



where 
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V(i) := ^Mm (12) 
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and is understood as a function of by inverting (fTUj). Specifically for the coupling 
function (jS)), the original and canonically transformed field are related in a simple manner 

-/acj) = sin(-\/o"0) for cr > and vH'/' — sinh(A/|cr|0) for cr < . (13) 

2.4 Transformation singularities 

Both the conformal and canonical transformation can fail to be one-to-one and thus become 
singular. For the conformal transformation this happens when d0/d0 = -F(0) = 0, i.e. at 

$c = = for < ^ < - 

V^(l - 60 6 

and when -F(0) diverges at 0c = Vv^ ^^ ^ > 0. While $c turns out to be a spurious 
singularity introduced by the transformation, 0c corresponds to a physical shear singularity 
at least when anisotropic perturbations are allowed jU]. 

The latter singularity clearly occurs for the canonical transformation, too, since now 
d0/d0 = /{(p)'^^"^ diverges at 0c. However, this is the only point where the transformation 
fails to be one-to-one and there is no analog of $c- The canonical transformation thus 
avoids the unphysical regime > $c of the conformal transformation. 



3 Singularities 

We can now investigate what this iniphes for singularities since we know the minimally 
coupled situation which is realized for the transformed variables. Then, for a regular 
potential the singularity corresponds to p = and ^ cx) jTU]. For a > 0, however, the 
potential can be singular, thus changing properties of space-time singularities [TT]. In the 
original variables, the behavior of the system can then be dramatically different for a < 
and 0" > as is evident from properties of the transformation. We will consider the former 
case first. 

3.1 Negative coupling constant 

For 0" < 0, the potential V{(j)) is regular everywhere such that in the transformed variables 
singularities occur in the usual form, p = and ^ oo. For purposes of quantization it 
will later be important to know what this implies for the original canonical variables p and 
since these are the variables that will be quantized. It follows from the preceding consid- 
erations that the behavior of will not be different from that of since the transformation 
(|T3|l is one-to-one. In fact, in the vicinity of the singularity we have oc exp(v^0), i.e. 
the original field also blows up. Interestingly, had the coupling function been of the form 
/(0) = 1 — cr0^^'^ for some positive e, however small, the integral in ()1()|1 would have been 
converging for ^ oo, that is a finite value of would have corresponded to a diverging 
and the possibility of a different type of singularity in the canonical variables would not 
have been ruled out. The case of a quadratic coupling function, which is preferred by phys- 
ical arguments since it does not require additional dimensionful constants, thus appears as 
the limiting case where singularities in non-minimal models occur as in minimal ones. 

It will also be of interest to see the behavior in the geometrical variables p in addition 
to 0. Using the relation between the scale factors (jHl), we conclude that p ^ implies 
p — >■ since /(0) 7^ 0. Thus, a singularity in the canonical variables implies a singularity 
in the geometrical variables occurring at the same values of the field and scale factor. The 
reverse conclusion is not true in general since p = together with ^ 00 could result in a 
finite value of p. However, since we already know that the canonical variables are regular 
except at p = 0, such cases are ruled out for a < 0. 

We are now ready to investigate the equations of motion generated by the transformed 
Hamiltonian (fTT|l . It turns out that in the vicinity of the singularity they can be solved 
exactly for many potentials (see also ^Hj for a general discussion in the minimally coupled 
case) : 

.~C [2k /., , „.„-./ 7x 7 TT ^ ~2i 



p = 2p- = ^ — ^Jn^ + 2pW{<j)) , = ^ , 7f = -pV,^(0) (14) 

where the dot denotes a derivative with respect to the conformal time 77 since we chose 
A^ oc a/Pi . In the first equation, we have eliminated the connection c, using the Friedmann 
constraint H ^ 0. We should now specify the form of the original field potential. The most 
common choices are the quadratic and quartic functions U{(f>) = \w?(t)'^ or f/(0) = tA0^. 



The corresponding transformed potentials then are 

^^ rr^tanhU^ and Vi4>) = A^^Mv^4>r ■ (15) 

^^' 2a cosh(v/^0)2 ^^' 4:a^ ^^ ^' ^ ' 

Note that these expressions are exact and work for any values of the field, which is one 
example for simplifications of the canonical compared to a conformal transformation. In 
the weak-field limit (a^^ <^ 1) one recovers the original potentials. On the other hand, 
when approaching the singularity (p -^ 0, acfP' -^ oo), the first potential vanishes, whereas 
the second one asymptotically tends to a constant, jXa"^. In both cases, the potential term 
in the p-equation (fT^ can be neglected as p ^ 0, while the righthand side of the yf-equation 
goes to zero. The latter implies that the field momentum is a constant, vr = ttq = const, 
which will thus dominate over p^V{(j)) sufficiently close to p = 0. We are, therefore, left 

with just two simple equations p = A/2K/37ro and (f) = ttq/p. Assuming that p{ri = 0) = 0, 
we get the solutions 

Piv) = T^oJ Y^ °^ ^ ' 4>{v) = J^^'^V^'^'^V- (16) 

Remarkably, the solutions do not distinguish between the two potentials and are recog- 
nized as those describing the massless scalar field. Note, also, that the asymptotic time 
dependence of the field does not depend on the initial field momentum. Converting back 
to the original variables, we write 



rs^ 



isinh(v/^0) = "^ " "L oc - ( - ) (17) 



P A /2/« V ^+^/^ 

We conclude the subsection, noting that in the limit of minimal coupling a,^ ^ (more 
precisely, a should go to zero faster than 0^^ to conform with -^ cxo) solutions ()17|) 
reduce to (fT^ . (In the first equation (fT7|) . one has to use the I'Hospital rule when taking 
(J to zero. For doing so, one must not ignore the first term with a positive power of t] as 
we did in the final step of this equation above at fixed a.) 

3.2 Positive coupling constant 

In this case, the mapping from to is not one-to-one. Compared to a conformal trans- 
formation, one such singularity at <l>c is removed which thus turns out to be spurious and 
just introduced by the choice of transformation. The remaining singularity at (pc, however, 
is different because the coupling function vanishes at this value, which thus corresponds to 
a singular point even in the original system. Nevertheless, it also appears as a point where 
the canonical transformation fails to be one-to-one. 



As a consequence, the effective potential V{(j)) diverges at this point. There is thus a 
potential wall which could simply reffect back when the critical value is approached and 
thus restrict the allowed values for to two disconnected regions. Whether or not this may 
happen in isotropic models, an embedding in anisotropic ones shows that 0c generically 
corresponds to a shear singularity rather than just a turning point of jH]. 

The irregular potential obtained after the canonical transformation thus implies a new 
type of singularity which occurs at a finite value of where /(0) =0. As follows from 
the considerations in JOj, this happens at a diverging and possibly non-zero value of p. 
This is a new type of singularity and presents an interesting test to singularity removal 
schemes of quantum gravity: Such schemes are usually based on quantum modifications 
at small scales and thus naturally arise around p = 0. But if variables are allowed to be 
of intermediate value at a classical singularity, it would be difficult for quantum gravity to 
remove such a singularity and at the same time preserve classical behavior on large scales. 

For a canonical quantization such as loop quantum gravity one uses canonical variables 
and it is their appearance at classical singularities which is relevant. Thus, we have to 
determine the values of p and at the new type of singularity when = 0c. While 
is finite there, the fact that p does not diverge together with /(0c) = implies that the 
canonical variable p must be zero, unlike the geometrical variable p. Thus, singularities 
always occur in regimes of small p where quantum gravity can dramatically change the 
behavior and has the potential to remove singularities. 

4 Loop quantum cosmology 

Loop quantum gravity is based on Ashtekar variables, defined by a densitized triad and a 
connection, which allow one to set up a kinematical quantum framework in a background 
independent manner. For minimal coupling, the densitized triad and connection are con- 
jugate, while for non-minimal coupling the densitized triad has to be multiplied by the 
coupling function in order to retain canonical variables. Upon reducing to isotropic met- 
rics, these canonical fields directly give the components p and c defined before, in addition 
to the scalar and its momentum vr. These are the variables that a loop quantization of 
the model should be based on since any canonical transformation such as that to and n 
may change properties and may not be representable by a unitary transformation at the 
quantum level. 

As in models with minimal coupling [7| ^1 E] , loop quantum cosmology in the triad 
representation is based on wave functions ■{/'(/i, 0) where /i is a label of eigenvalues of p. 
The classical Friedmann equation for p is then replaced by a difference equation for the 
wave function of the type 

Ci{jl)i){fx + 1, 0) - 2Co{jx)ll){jl, 0) + C-i{jl)i){fl - 1, 0) OC i^mattcr(/i)V'(/i, 0) (18) 

with the matter Hamiltonian operator i/mattcr(/i) which for standard matter is diagonal 
on states |/i). In the limit of large scale factors /i ^ 1, this difference equation yields 
the Wheeler-DeWitt differential equation while at small scales the discreteness of quantum 
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geometry is essential. The criterion for non-singular behavior (see also ^1] for details) then 
consists in unique extendability of the wave function even across classical singularities for 
specific values of ft. In the standard case, singularities are reached for fi = 0, where 
the kinematical structure provides a new branch since there are two sides to the classical 
singularity corresponding to opposite orientations of the triad. While they cannot be 
connected classically, it is possible to extend the wave function uniquely. To see this, one 
has to use a more detailed form of coefficients of the difference equation and discuss when 
they become zero. It turns out that, although leading coefficients C±i can become zero, 
this does not spoil the recurrence scheme of the difference equation. The wave function 
is then uniquely extended and quantum gravity can tell us about the other side beyond 
singularities ^E^. 

For this argument it is important to provide a candidate for the other side, such as /i < 0, 
and then show that dynamics extends the wave function uniquely. Thus, this scheme only 
works when classical singularities are indeed located at values of vanishing fi (or degenerate 
triads) since there would be no new region otherwise. While this is automatically the 
case in isotropic models with standard types of matter fields, it is non-trivially realized in 
anisotropic models JHl or even spherical symmetry ^7] . In anisotropic models, for instance, 
the Kasner behavior where one metric component always diverges at a singularity suggests 
a general behavior where not all metric components vanish at singularities. However, it 
turns out that densitized triad components which underlie the loop quantization always 
vanish at anisotropic classical singularities. Thus, the same argument to remove classical 
singularities applies. This is a general scheme, but non-minimally coupled models were 
not included so far. In fact, a second crucial ingredient is the assumption of a matter 
Hamiltonian operator diagonal on triad eigenstates. Classically, this corresponds to a 
matter Hamiltonian which depends only on spatial geometry but not on curvature. If 
there are curvature terms and -^matter is not diagonal in geometric variables, the right 
hand side of ()18|) could also contribute terms proportional to ?/'(/i ± 1, 0) which would 
change the recurrence and introduce new places where the leading coefficients vanish. The 
recurrence would then have to be re-analyzed and singularity removal would be less general 
since the coefficients would depend on the matter field through the Hamiltonian, and not 
just on ft. 

We can now test this scheme with what we have learned about non-minimally coupled 
isotropic models. First, there may be singularities where even the geometric triad variable 
p may not vanish. However, it turned out that the canonical variable p does vanish. 
Since the difference equation is based on the quantization of this variable, non-minimally 
coupled models are included in the non- singularity scheme of loop quantum cosmology. As 
the discussion here illustrates, this happens in a way depending on dynamical details of 
the models. 

At first sight, one could also expect that curvature couplings of a non- minimal model 
lead to non-diagonal matter Hamiltonians since there are curvature components in the 
matter terms. Here, the previous expressions show that the matter Hamiltonian © in 
canonical variables does depend only on p in addition to the matter fields, but not on the 
connection component c which determines curvature. Thus, the matter Hamiltonian even 



with non-miniinal coupling is diagonal in the canonical triad representation and the general 
singularity removal scheme applies. This remains true for positive spatial curvature not 
spelled out explicitly here, since we would only have to replace ?? by c^ + 1 [iH]. Thus, 
non-minimally coupled isotropic models are non-singular in loop quantum cosmology in 
the same manner as minimally coupled ones. 

5 Conclusions 

Applying a canonical transformation to a non-minimally coupled isotropic cosmological 
model we have shown that for classically allowed coupling functions singularities occur as 
in the minimally coupled case. For quadratic and quartic scalar potentials, moreover, the 
asymptotic approach to a singularity is as in the massless, free minimal case. 

This has direct implications for the question of singularity removal in quantum gravity. 
Loop quantum cosmology [13j has provided a scheme of singularity removal which has 
been shown previously to apply to isotropic [inj, diagonal class A Bianchi models [T?)|IT!l]. 
and spherically symmetric and polarized cylindrical wave models |17| 1^ . In all these 
cases, it was shown that quantum dynamics provides a well-posed recurrence scheme for 
wave functions across classical singularities. It was therefore important to know at which 
geometrical configurations a classical singularity occurs, such as at vanishing scale factor 
p = for isotropic models coupled to a scalar. All these arguments were independent of 
the form of the matter Hamiltonian and thus present a quantum geometry effect, provided 
that the classical matter Hamiltonian did not depend on curvature. 

For the purpose of quantum gravity, curvature couplings fall into two different classes, 
depending on whether they require a non-zero Planck length or not. Most curvature 
coupling terms in effective matter actions require dimensionful coefficients with a new 
length scale to compensate the dimensions of curvature factors. Such terms cannot arise 
in the classical Hamiltonian used to set up a quantization. Such actions with curvature 
coupling terms will thus not be loop-quantized; they would rather follow from a quantum 
theory of gravity such as loop quantum gravity in an effective description, see e.g. pT] . 
The curvature coupling of a non-minimal scalar, on the other hand, is allowed classically 
provided that the coupling function is quadratic in 0, but it was not covered by previous 
discussions of singularities in the context of loop quantum gravity. 

With the results of the present paper, we can fill this gap. There are two potential 
problems: the Hamiltonian could be curvature dependent, and different types of singular- 
ities could occur. For the classically allowed form of coupling function which is quadratic, 
we have seen that classical singularities always occur at vanishing p = f{4>)p which is one 
of the canonical variables, in the same form as in minimal models. Moreover, the form 
of quantum dynamics near classical singularities is the same whether a scalar is coupled 
minimally or non-minimally. The argument of singularity removal is thus unchanged. 

The discussion of this paper shows that this conclusion is realized quite non-trivially: 
The classically allowed case of a quadratic coupling function is just the limiting case where 
an infinite scalar is mapped to an infinite by the canonical transformation. For coupling 
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functions of higher degree, the possibihty of singularities at finite values of and non-zero 
p cannot easily be ruled out. In such a case, the usual singularity removal mechanism and 
the resulting non-singular picture of quantum geometry would have been very different, if 
applicable at all. Since higher degrees than two in the coupling function can only arise 
effectively, but not in a classical action yet to be quantized, the singularity structure rele- 
vant for isotropic loop quantum cosmology does not change for non-minimal coupling. As 
with calculations of black hole entropy [22] we have seen that, although the non- minimal 
situation appears initially crucially different from the minimal one, the methods of loop 
quantum gravity are general enough to encompass automatically also non-minimal situa- 
tions in quite non-trivial ways. 
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